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Abstract 

A detailed general relativistic (GR) calculation of the neutrino path for a general metric de- 
scribing a rotating star is studied. We have calculated the neutrino path along a plane, with 
the consideration that the neutrino does not at any time leave the plane. The expression for the 
minimum photosphere radius (MPR) is obtained and matched with the Schwarzschild limit. The 
MPR is calculated for the stars with two different equations of state (EOS) each rotating with two 
different velocities. The results shows that the MPR for the hadronic star is much greater than 
the quark star and the MPR increases as the rotational velocity of the star decreases. The MPR 
along the polar plane is larger than that along the equatorial plane. 
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I. INTRODUCTION 



Gamma Ray Bursts (GRB), the possible engines for GRBs and its connection with the 
neutrino production is a field of high current interest. It was proposed that the neutrino- 
antineutrino annihilation to electron-positron pairs in compact stars is a possible and im- 
portant candidate to explain the energy source of GRBs. The previous calculations of the 
reaction vv — ■* e + e~ in the vicinity of a neutron star have been based on Newtonian gravity 
[l[ [2], i.e. (2GM/c 2 R) << 1. The effect of gravity was incorporated in refs. 0,0], but only 
for a static star. 

Neutron stars are objects formed in the aftermath of supernova. The central density of 
these stars can be as high as 10 times that of normal nuclear matter. At such high density, 
any small perturbation, e.g. spin down of the star, may trigger the phase transition from 
nuclear to quark matter system. As a result, the neutron star may fully convert to a quark 
star or a hybrid star with a quark core. It has been shown [5j that such a phase transition 6] 
produces a large amount of high energy neutrinos. These neutrinos (and antineutrinos) could 
annihilate and give rise to electron-positron pairs through the reaction uu — > e + e~ . These 
e + e~ pairs may further give rise to gamma rays which may provide a possible explanation 
of the observed GRB. Furthermore, the rotating neutron star has been shown [?] to produce 
the observed beaming effect. At present, it is necessary to have a better understanding 
of the energy deposition in the neutrino annihilation to e + e~ in the realistic neutron star 
environment. 

We would like to study the v + V — ¥ e + + e~ energy deposition rate near a rotating 
compact star. This reaction is important for the study of gamma ray bursts. The General 
Relativistic (GR) effect increase the efficiency of the process immensely, but the inclusion 
of the rotational effect is yet to be incorporated and studied. The geodesic of neutrinos 
are also important in the study of pulse shapes and accretion disc illumination [8] to name 
a few. Therefore the path of neutrino (or generally of massless particle) is of immense 



importance and needs a detailed study. The neutrino path for Schwarzschi 
the equatorial pla ne can be found in text books [9] and different papers 



Fukuyama 
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d metric along 
3|. Asano and 



ll| did the same calculation near a thin accretion disc using Kerr metric. 



Prasanna and Shrubabati 



121 ] studied it for a slowly rotating star. To address all the above 



problems we present a detailed GR study of the neutrino path for a most general metric 
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describing a rotating star [13j along a plane. We have made our calculation using two 
different EOS, one quark and the other hadronic. 

In this paper first we will discuss about the metric, the EOS and the star structure. Next 
we will present the detailed GR calculation of the neutrino path and minimum photosphere. 
Finally we will present our results for the two EOS and have a brief discussion. 

II. THE STAR 

The structure of the star is described by Cook-Shapiro- Teukolsky (CST) metric [3] 

ds 2 = -e t+p dt 2 + e 2a {dr 2 + r 2 d6 2 ) + e~<- p r 2 sin 2 6{d<j) - udt) 2 . (1) 
Accurate models of rotating neutron stars for tabulated EOS can be computed numerically 



using the 'rns' code 
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161 ] . This computer code computes the metric functions a,j,p 
and io appearing in the axisymmetric metric, and these metric functions depends only on the 
coordinates 9 and r. The metric function u is the term responsible for the frame dragging 
effect and would vanish if the rotational velocity (fl) is zero. The coordinate r is related to 
;he standard radial coordinate that appears in the Schwarzschild metric, r s , by r s = re^~ p ^ 2 



171 ] . In the limit of zero rotation, the following combination of metric functions are 



Lim^ re p = (2) 



I _ 2M 

r s 

Um^e^ 2 = yjl - ^ (3) 
Um^e-^^dr = I ^ E . (4) 

We have previously mentioned that tabulated EOS are needed to compute the code 
numerically. In this paper we have computed for two different EOS, the quark EOS and the 
ladronic EOS. The hadronic EOS has been evaluated using the nonlinear Walecka model 
1^ | . The Lagrangian density in this model is given by: 
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i 

- jP^Pf + \ m2 p PlPa - \bm N ( g(7N af - -c( grTN a) 4 

+ ^{i'fdn - m e )il> e (5) 

The Lagrangian in eqn. 5 includes nucleons (neutrons and protons), electrons, isoscalar 
scalar, isoscalar vector and isovector vector mesons denoted by tpi, ip e , a, uj 11 and p a,fl , re- 
spectively. The Lagrangian also includes cubic and quartic self interaction terms of the 
a field. The parameters of the nonlinear Walecka model are meson-baryon coupling con- 
stants, meson masses and the coefficient of the cubic and quartic self interaction of the a 
mesons (b and c, respectively). The meson fields interact with the baryons through linear 
coupling. The u and p meson masses have been chosen to be their physical masses. The 
rest of the parameters, namely, nucleon- meson coupling constants (^-, and ^-) and the 
coefficients of cubic and quartic terms of the a meson self interaction (b and c, respectively) 
are determined by fitting the nuclear matter saturation properties, namely, the binding en- 
ergy/nucleon (-16 MeV), baryon density (po=0.17 fm~ 3 ), symmetry energy coefficient (32.5 
MeV), Landau mass (0.83 m n ) and nuclear matter incompressibility (300 MeV). We have 
used a stable three-flavour quark matter EOS, obtained from the standard Bag model with 
B 1 / 4 = UbMeV. 

The shape of a fast rotating neutron star becomes oblate spheroid fl3 |. The star gets 
compressed along the z-axis and along x and y-axes, it bulges by equal amounts, the polar 
radius is thus smaller than equatorial radius. 

III. GR CALCULATION 

The metric is independent of 't' and '4>', the coordinates are cyclic, hence the correspond- 
ing covariant generalized momenta is constant [9j , i. e 

Pt = Po = const. = —E 

P<t> — Ps — const = L (6) 
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The magnitude of the 4 vector energy momentum is given by 9|, 

gijpY + /i 2 = o, 

where \i is the rest mass of the particle. Writing it explicitly we have 

fi'ooP 02 + 9\\V l2 + 922P 22 + 933P 32 + 93oP 3 P° + 9o3P°P 3 + /i 2 = 0. 



(7) 



(8) 



These are contravariant momenta, and to find the contravariant momenta p\ we calculate 
the inverse matrix , of the metric given in eqn. 1 



_ e -(7+p) 




-ue^ +p) 





e" 







-2a 








\e~ 2a 



-use 



-(rr+p) 










-(u 2 e~^ - ^~ P) 



r' 2 sin 2 6 



Using the original g^ u and the inverse g pv matrix we calculate the contravariant momenta, 
which are given by 



P 



9°°Po + 9° 3 P3 = e~^ +p \E - uL) = B(co) 



P 3 = 9 3 °Po + 9 33 P3 



e-^ +p) \uB(uj) + 



,2p 



r 2 sin 2 ( 



(9) 



We consider the particle motion is in a particular 9 = const plane, and orient the coor- 
dinate system such that the particle lies in the equatorial plane for 9 — |- . The particle has 
at start, and continues to have zero momentum in the given plane i.e p e = p 2 = 0. 
Finally, substituting these values in the above eqn. 8 we get 

L 2 



e -h+p) B 2 



^™e^ + e 2 «& 2 + // = (). 
r z svn z 9 dX 



(10) 



if, A = -i.e propertime per unit rest mass 
p 



dr dr dr 



dr 



dX dr dX ^ dr 
dr ^dr d(p^ ^dr ^d<f)^dX^ 

dX dcf) dr d<p dX dr 

dX y dd> )P 



(11) 



We define E — — and L — -. As the particle is massless (neutrino), therefore we define 



where b is the impact parameter, substituting this in eqn. 10, we have 



e 2 «A>(l-u;&) + ^^] 2 
d(p r 2 sin z 



(l-ujb) 2 + 



b 2 



r 2 sin 2 9 



; 7+3p = Q> 



(13) 



The lagrangian of the system we are considering is given by 
where, x 1 — A the affine parameter. Therefore writing it explicitly, we get 



(14) 



L = (-e 7+p + e 1 ~ p uj 2 r 2 sin 2 6)i 2 + e 2a r 2 + e^- p r 2 sin 2 6^ 2 - 2e~ 1 ~ p ujr 2 sin 2 6&. (15) 

Using the above lagrangian we can write the covariant momentums, and they are given in 
terms of total energy and total angular momentum 

dL 



p = pt = = {- e 1+p + e^ p uj 2 r 2 sin 2 9)i - e"<~ p ujr 2 sin 2 6<f) = -E 
at 



P3 =P<P 



dL 



e 7 p r 2 sin 2 6(j) — e 1 p ur 2 sin 2 9t = L. 



(16) 



Having written the momenta in terms of total energy and total angular momentum it is 
quite simple to solve for t and <fi. They are given by, 

E-uoL 







+ 



t = 

uj{E — ujL) 



(17) 



ei-Pr 2 sin 2 9 e^+P 
The angle 6 r between the particle trajectory and the tangent vector to the orbit can be 
derived by constructing the local lorentz tetrad k p v for our metric 



e (7+P)/2 




—e^' p ^ 2 ujrsin8 
The angle 6 r is given by 





e a 





U 1 

tanOr = T77T 







re 







e^-p)/ 2 rsin9 



k\V r 



v 2 kW* + k^v* 

where, V r = j and — | are the local velocities. Using eqn. 17 can be written as 



(18) 



y<t> 



,;2p 



E — ujL r 2 sin 2 9 



+ u = A(r,6) +uj. 



(19) 



Therefore the angle 9 r is 



e a A(r,9) +uo dr 



e (i-p)/2 rs i n Q' A(r,\ 
Squaring the above equation and writing as 



(20) 



>- - 1 -tan 2 6 r . (21) 



2a 



we get the final form of Substituting this value in eqn. 13, we get 

2a , A(r,6) 2 e^- p ^r 2 sin 2 e 2q e 2p b 



tan 6 r [uj(l — cob) + 



L A(r,6) +U 1 e 2a L v ' r 2 sin 2 9 i 

-e^ +p \l - uob) 2 + e ^+ 3 ^ = 0. (22) 

r 2 sin 2 9 

This equation can be solved using the potentials obtained from rns code to obtain a min- 
imum radius r = R, the minimum photosphere radius, below which a massless particle 
(neutrino) emitted tangentially to the stellar surface {9r = 0) would be gravitationally 
bound. 

In the limit in which CST metric reduces to Schwarzschild metric, i.e 



Lim^ore p = 

eqn. 22 (for 9 = 7r/2) reduces to the equation for b obtained by Salmonson and Wilson 
i.e 

r q cos9 r 



a, 



l _ 2M 



IV. RESULTS 

The minimum photosphere is calculated solving eqn. 22 using the potential functions 
obtained from the 'rns' code. Starting our calculation by choosing the central energy 
density of the star to be 1 x 10 15 gm/cm 3 . Fig. 1, 2, 3, 4 and 5 gives the nature of the 
equations of state used. In the figures x — cos9. Fig. 1 shows that the quark matter 
EOS, considered here is much steeper than the hadronic matter EOS. Fig. 2 (for quark 
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energy density (gm/cm ) 



FIG. 1: Variation of pressure with energy density for quark and hadronic matter EOS. 

matter EOS) shows that at the centre of the star the pressure is maximum, and as we go 
outside it falls off, and becomes zero outside the star. Along the pole the pressure falls off 
in a much steeper way than along the equator, as it is of much shorter length. As the the 
rotational velocity decreases the equatorial radius of the star decreases but the polar radius 
increases (although still less than equatorial radius). For the keplerian velocity the star is 
maximally deformed and as the rotational velocity of the star decreases the star regains a 
much spherical shape. Fig. 3 show the same nature for a hadronic star. In fig. 4 and 5 the 
variation of energy density is shown for the quark and hadronic star simultaneously, and its 
nature is more or less similar to that of the behaviour of pressure discussed above. 
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FIG. 2: Variation of pressure along the radial direction of the star for two different rotational 
velocities each with two different values of x f° r the quark matter EOS. 
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FIG. 3: Variation of pressure along the radial direction of the star for two different rotational 
velocities each with two different values of x f° r the hadronic matter EOS. 
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FIG. 4: Variation of energy density along the radial direction of the star for two different rotational 
velocities each with two different values of x f° r the quark matter EOS. 




FIG. 5: Variation of energy density along the radial direction of the star for two different rotational 
velocities each with two different values of x f° r the hadronic matter EOS. 
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Using the quark matter EOS on the 'rns' code, the keplerian velocity of the quark star 
comes out to be 0.89 x 10 4 s _1 . For comparison we have also computed the 'rns' code with 
rotational velocity of 0.5 x 10 4 s _1 . The same treatment is done also for the hadronic EOS 
where the keplerian velocity is 0.61 x 10 4 s~ 1 and for comparison the other rotational velocity 
was chosen to be 0.4 x 10 4 s _1 . The code solves the metric for the given EOS and gives the 
different potential functions as a function of r and 9. Solving eqn. 22 with these values of 
potential functions for different 9 we obtain the value of minimum photosphere for different 
planes. Table 1. sums up all our results in a compact form. 

Let us now analyze the table given above. It points out the fact that as the rotational 
velocity decreases the mass of the star also decreases. For the same central energy density, 
the mass of the quark star is much greater than that of the hadronic star but the radius 
is much smaller. It signifies that the quark matter EOS, considered in our work is much 
steeper than that of hadronic matter EOS as pointed out previously in the figures. As the 
rotational velocity decreases the equatorial radius decreases but the polar radius increases. 
It shows that the star is maximally deformed for the keplerian velocity and as the rotational 
velocity of the star decreases it tries to regain a more spherical shape. A static star is of 
spherical shape, where polar and equatorial radius are same. 
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In the above table we have tabulated the minimum photosphere radius (MPR) for three 
different values of %■> *.e for three planes. Along the equator (x = 0), along the pole 
(x = 0.99) and along a plane lying at \ — 0-5. The MPR is minimum along the equator 
and maximum along the pole. For \ = 0.5 it lies somewhere in between these two values. 
The MPR is minimum for the quark star rotating with keplerian velocity and is maximum 
for the hadronic star rotating with 0.4 x 10 4 s _1 velocity. The MPR is much greater for 
the hadronic star than the quark star. As the rotational velocity decreases the MPR shifts 
outward from the centre of the star toward the surface. 



V. SUMMARY AND DISCUSSION 

In this paper we have addressed the problem of path of the neutrino and the radius of 
minimum photosphere. We have done a complete GR calculation of the neutrino path for 
the most general metric describing a rotating star, and have obtained its geodesic equation 
along a given plane. We have calculated the MPR for four cases, i.e stars with two different 
EOS and both rotating with two different velocities. Previous calculation of the neutrino 
path was either done for a static star 3|, |4| or for a slowly rotating star [12J and only along the 
equatorial plane. We have shown that our results also matches very well with the previous 
findings [3( for the Schwarzschild limit. We have found that the MPR is maximum along 
the pole and minimum along the equator. The MPR is much greater for the hadronic star 
than that of the quark star. As the rotational velocity decreases the MPR increases and 



is maximum for the static star. Prasanna and Shrubabati [12] had showed that the MPR 
is inversely proportional to the rotational velocity of the star. Salmonson and Wilson 
had shown that for the static star the MPR limit is R = 3M, and that is very close to the 
surface. So our results are at par with the previous findings in those limits and goes beyond 
them. 

Finally we would like to mention that this calculation of neutrino path is very important in 
the sense that this forms the heart of different problem like GRB central engine, pulse shape 
and accretion disc illumination. This path is the general path followed by any massless 
particle (photon) in the vicinity of a compact object. So the calculation might also be 
important to those problems. Currently we are trying to address other problems related to 
neutrino path. 
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